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INTRODUCTION 

The  overall  objective  of  the  research  described  In  this  final  report  is 
to  develop  diagnostics  and  models  required  for  measurement  and  prediction  of 
the  properties  of  electrical  discharges  in  molecular  gases.  Such  discharges 
are  important  to  the  Air  Force  because  of  their  use  in  devices  such  as  high 
power  switches,  negative  ion  sources  and  plasma  processors  and  because  of 
their  role  in  phenomena  such  as  lightning  and  corona.  The  research  program 
includes;  a)  the  development  and  application  of  laser  fluorescence  and 
absorption  techniques  to  the  measurement  of  excited  state  densities  in 
electric  discharges,  b)  the  measurement  of  the  properties  of  the  excited 
states  of  importance  in  molecular  discharges  of  current  interest,  c)  the 
application  of  various  diagnostic  techniques  to  the  determination  of  the 
characteristics  of  moderate  current  density,  transient  electrical  discharges, 
and  d)  the  quantitative  comparison  of  these  results  with  appropriate 
discharge  models. 

During  the  past  year  we  have  (a)  completed  measurements  of  the 

collisional  destruction  of  metastable  molecules  in  various  rotational  and 

vibrational  states,  (b)  initiated  measurements  of  excitation  coefficients  for 

Hz ( c  II u)  metstable  molecules  by  discharge  electrons,  and  c)  attempted  to 

model  experimental  determinations  of  the  time  dependence  of  00(a^A^) 

metastable  densities  in  pulsed  electrical  discharges.  The  properties  of 

metastable  Hz  are  of  current  interest  because  of  the  use  of  discharges  in  Hz 

as  switching  devices  ( thyratrons ) ,  as  negative  ion  sources,  and  as  infrared 

lasers.  The  behavior  of  0o(*A  )  metastables  in  electric  discharges  is  of 

1  g 

interest  because  of  the  possibility  of  discharge  lasers  utilizing  excitation 


-1- 


transfer  from  the  metastable  to  I  atoms  which  can  undergo  stimulated 
emission  in  the  near  infrared. 

3 

The  results  of  work  on  the  c  nu  state  of  H2  is  summarized  in  Sec.  II. 
The  diagnostic  techniques  and  some  initial  results  connected  with  the 
measurement  of  electron  excitation  coefficients  for  the  ^(c^n  )  state  are 
discussed  in  Sec.  III.  The  results  obtained  thus  far  in  the  analysis  of  the 
transient  O2  densities  in  pulsed  discharges  are  discussed  in  Sec.  IV. 
Recommendations  for  future  work  are  made  in  Sec.  V. 


SECTION  II 


H2  METASTABLE  DESTRUCTION 

In  this  section  we  summarize  our  measurements  of  rate  coefficients  for 

3 

the  destruction  of  various  rotational  and  vibrational  levels  of  the  c  nu 
metastable  state  of  H2  in  collisions  with  H2.  We  also  summarize  a  brief 
search  for  collisional  coupling  between  H2(c  II u)  levels  and  nearby  aJ£^ 
levels.  The  experimental  technique  and  the  results  for  the  N  =  1  rotational 
level  the  v  =  2  vibrational  state  have  been  published.*  We  will  therefore 

only  briefly  review  the  experimental  technique  before  presenting  the  results. 

3 

The  lower  rotational  and  vibrational  levels  of  the  c  Jlu  state  and  of  the 

nearby  aJEg  state  of  H2  are  shown  in  Fig.  1.  The  levels  investigated  in  this 

work  are  marked  with  an  asterisk.  As  described  in  Ref.  1,  the  relative 

excited  state  densities  were  measured  by  observing  the  absorption  of  radiation 

in  transitions  to  the  g^£  +f  i^jj  and  j"*A  states.  Figure  2  shows  the 

8  o  & 

experimental  arrangement  and  Fig.  3  shows  a  representative  absorption 

transient.  As  in  our  previous  experiments,  the  metastable  density  was 

suddenly  reduced  by  exciting  a  large  fraction  of  the  metastables  to  a  higher 

state  from  wnich  a  significant  fraction  could  radiate  to  the  b-*l  +  state  and 

u 

then  dissociate.  The  time  constant  of  the  recovery  of  metastable  population 
was  determined  by  fitting  an  exponential  to  the  absorption  waveform  during  the 
recovery  period  as  shown  in  the  expanded  portion  of  the  absorption  transient 
shown  in  the  lower  part  of  Fig.  3.  The  principal  changes  in  the  experimental 
apparatus  are  the  use  of  a  balanced  detection  system  so  as  to  reduce  the 
effects  of  intensity  fluctuations  in  the  output  of  the  cw  dve  laser  source 
used  for  the  absorption  measurements  and  the  replacement  of  the  N2~punp  laser 


Nd:YAG  pumped 
pulsed  dye  laser 


trig. 


H.V. 

pulser 


SIGNAL  <mV> 


Figure 


Upper  trace:  overall  transient  absorption  waveform  showing 

metastable  increase  during  discharge  pulse  and  depletion  caused  by 

pulsed  laser.  Lower  trace:  expanded  absorption  waveform  showing 
3 

recovery  of  c  nu  (N“l,  v*0)  metastable  population  and  exponential 


fit  to  data. 


for  the  pulsed  dye  laser  with  a  much  more  powerful  and  reliable  Nd:YAG  pump 


laser. 


Representative  rate  coefficient  data  for  the  N  =  1,  v  =  0  level  are 
plotted  in  Fig.  4  as  a  function  of  current  for  various  discharge  pressures. 

The  straight  lines  are  a  least  squares  fit  to  the  data  assuming  a  linear 
dependence  of  the  metastable  destruction  on  the  discharge  current  and  on  the 
H2  density.  The  rate  coefficients  for  H2(c^nu)  destruction  at  zero  discharge 
current  are  plotted  as  a  function  of  vibrational  quantum  number  v  for  N  =  1  in 
Fig.  5  and  as  a  function  of  rotational  quantum  number  for  v  =  1  in  Fig.  6. 

The  error  bars  are  indicative  of  the  uncertainty  in  the  average  values  at  a 
given  current  and  pressure.  The  results  of  these  experiments  are  also 
tabulated  in  Table  I.  The  present  results  for  the  v  *  2,  N  **  1  level  are  in 
very  good  agreement  with  our  published  results.*  We  note  that  there  is  very 
little,  if  any,  change  in  the  destruction  rate  coefficient  with  vibrational 
quantum  number  and  only  about  a  102  change  with  rotational  quantum  number.  In 
view  of  the  very  large  magnitude  of  these  quenching  rate  coefficients  and 
associated  cross  section  (~80  A^),  each  gas  kinetic  collision  results  in 
quenching  and  we  do  not  expect  much  change  with  initial  level. 


A  brief  effort  was  made  to  determine  whether  the  collisional  quenching  of 
the  c3nu  levels  resulted  in  excitation  transfer  collisions  with  nearby 
levels.  In  these  experiments  the  population  of  one  of  the  more  highly 
absorbing  levels  of  the  c  nu  state  was  perturbed  by  the  pulsed  laser  and  the 


*1  + 

resulting  change  in  the  a  level  density  was  observed  by  tuning  the  cw 

laser  to  the  appropriate  a3*;  +  to  d3n  transition.  Table  II  shows  the 

g  “u 

1  1  X 

combinations  of  c  Hu  and  a  levels  observed,  the  absorption  line  used,  and 

whether  or  not  the  a3£^+  population  was  observed  to  decrease.  Figure  7  shows 


the  observed  decreases  in  absorption  for  the  v  =  1,  N  =  1  level  and  for  the 


WWC WWW 


Figure  4.  Measured  and  calculated  reciprocal  decay  constants  for  v  =  0,  N  * 
1  levels  of  the  H2^c^u^  state  vs  discharge  current  for  various 


pressures. 


Figure  6.  Measur 


rate  coefficients  vs  rotational  quantum 


Table  I 


Measured  rate  coefficients  for  colllslonal 
destruction  of  vibrational  and  rotational 
levels  of  H2(c  nu)  metastables 


Vibrational 

level 

Rotational 

level 

Rate 

coefficient 

V 

N 

(10~ 

15 

3  - 

m  s 

0 

I 

1.85 

± 

0.11 

l 

l 

1.91 

± 

0.10 

1 

2 

1.92 

± 

0.08 

1 

3 

1.70 

t 

0.12 

2 

1 

1.98 

t 

0.10 

3 

l 

1.93 

t 

0.10 

Table  II 


Levels  examined  in  search  for  collisional  coupling 

3  + 


between  C^nu  and  a  states  of  H2 


N 

V 

X  (nm) 

N 

V 

X  (nm) 

Signal 

observed 

c%u  (depleted) 

a3Ig+  (observed) 

1 

1 

588.82 

1 

1 

608.08 

* 

yes 

1 

0 

593.14 

1 

0 

601.83 

* 

yes 

1 

0 

593.14 

1 

1 

608.08 

no 

2 

0 

596.35 

1 

0 

601.83 

no 

c3n 

(depleted) 

3 

c  n  (observed) 

1 

1 

602.13 

3 

1 

588.39 

no 

The  observation  of  a  signal  does  not  necessarily  mean  that 
collisional  coupling  was  responsible.  See  text. 
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3  + 

v  *  0,  N  =  1  level  of  the  a  Eg  state  when  the  corresponding  levels  of  the 

3 

c  nu  state  were  depleted  with  the  pulsed  laser.  As  indicated  in  Table  II,  no 

other  signals  indicative  of  collisional  coupling  were  detected  that  were  above 

the  relatively  high  noise  level  resulting  from  the  short  radiative  lifetime 
3  +  2 

for  the  a  Eg  state.  We  have  not  completed  our  analysis  of  these 
waveforms.  Only  that  for  the  v  ■  1,  N  ■  1  level  has  the  expected  time 
dependence,  i.e.,  a  relatively  slow  decrease  of  a  state  population  and  a 
recovery  in  a  time  comparable  with  that  for  the  c^nu  levels.  Since  the 
waveform  for  the  v  =  0,  N  =  1  level  resembles  those  for  the  c  nu  levels, 
further  tests  for  proper  line  identification  are  necessary  before  final 
conclusions  can  be  reached.  Table  II  also  shows  that  for  the  v  *  1  levels  of 

3 

the  c  jiu  state  we  did  not  observe  a  decrease  in  the  density  of  the  N  =  3  level 
when  the  N  ■  1  level  population  was  depleted  with  the  pulsed  laser.  This 
observation  suggests  that  collisions  resulting  in  changes  of  N  by  2  are  not 
probable.  Changes  of  N  by  1  with  the  expected  absence  of  change  of  ortho-para 
symmetry  would  be  difficult  to  detect,  since  the  final  level  predissociates 
rapidly  (~10  ^  s).  These  observations  are  consistent  with  our  expectation 

3 

that  the  quenching  of  the  c  nu  metastables  is  too  rapid  to  be  accounted  for  by 
rotational  excitation. 


SECTION  III 


h2  metastable  excitation  by  electrons 

The  objective  of  these  measurements  is  to  determine  the  coefficients  for 

3 

the  excitation  of  H2(c  nu)  metastables  by  electrons  under  discharge 
conditions.  The  technique  being  developed  is  to  combine  measurements  of  the 
absolute  metastable  density  with  our  previous  measurements  of  the  collisional 
destruction  rate  coefficients  to  calculate  the  production  rate  and  excitation 
coefficients.  Comparison  of  these  experimental  excitation  coefficients  with 
theoretical  values  serves  to  test  proposed  excitation  cross  sections  used  in 
modeling  discharge  devices  such  as  positive  and  negative  ion  sources, 
thyratrons,  and  plasma  processors.  Since  the  work  described  in  this  section 
is  new  and  incomplete,  the  results  cited  illustrate  the  techniques  being 
developed  and  tested  and  should  not  be  regarded  as  final. 

The  apparatus  used  in  these  measurements  is  the  same  as  that  shown  in 
Fig.  2,  except  that  the  discharge  tube  was  modified  by  the  addition  of  two 
wall  probes  and  a  second  anode  for  the  differential  measurement  of  electric 
field  strengths  and  by  moving  the  cathode  out  of  the  optical  absorption 
path.  The  change  of  the  cathode  position  followed  the  discovery  that 
significant  absorption  in  the  discharge  tube  of  Fig.  2  occurred  in  the  cathode 
region.  A  photomultiplier  mounted  on  a  moveable  platform  was  added  to  allow 
observation  of  the  time  and  spatial  dependence  of  the  light  emission  while 
testing  for  discharge  nonuniformities  such  as  striations.  The  discharge  is 
operated  in  the  pulsed,  high  current  (3  ys  long  and  1  A  current)  mode,  since 
attempts  to  obtain  useful  absorption  without  striations  using  low  current  dc 
discharges  were  unsuccessful.  The  electric  field  probes  were  constructed  of 
0.08  mm  tungsten  wire  in  the  form  of  a  ring  adjacent  to  the  tube  wall  so  as  to 


minimize  perturbation  of  the  discharge.  High  voltage  oscilloscope  probes  was 
used  to  measure  the  time  dependent  floating  potentials  of  the  probes.  The 
impedance  of  the  probe  circuit  was  varied  to  demonstrate  lack  of  significant 
loading  of  the  plasma  by  the  probe.  Figure  8  shows  representative  discharge 
current,  probe  voltages,  and  absorption  transients.  We  see  that  the  discharge 
current  has  reached  a  plateau  by  about  I  ps  after  the  initiation  of  the 
discharge.  By  reading  the  current,  probe  voltage  differences,  and  absorption 
at  this  time  we  were  able  to  make  measurements  before  the  build  up  of 
significant  longitudinal  spatial  variations  in  the  discharge  light  output  and, 
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presumably,  before  the  build  up  of  longitudinal  electric  field  variations. 

The  solid  points  of  Fig.  9  show  the  resultant  E/n  vs  nR  data  for  these 
pulsed  discharges.  The  solid  curve  shows  the  predictions  of  theory^  when  the 
charged  particle  losses  are  assumed  to  be  determined  by  ambipolar  diffusion 
and  when  the  ionization  coefficients  are  calculated  using  the  cross  section  of 
Buckman  and  Phelps.  The  open  points  show  the  measured  E/n  vs  nR  values 
calculated  from  the  pulsed  discharge  data  of  Breusova^  and  the  dc  positive 
column  data  of  Giintherschulze. ^  On  the  basis  of  the  current  dependence  found 
by  Breusova  we  expect  the  steady  state  values  of  Giintherschul ze  to  be  low. 

Note  that  our  discharge  pulse  lengths  are  much  shorter  than  those  of  Breusova 
so  that  we  can  operate  at  higher  currents  than  those  at  which  she  found 
significant  heating  effects.  Although  the  agreement  in  magnitude  between  our 
experimental  results  and  theory  is  good,  the  variation  of  the  experimental  E/n 
with  nR  is  too  slow  at  the  larger  nR  values.  This  divergence  may  he 
associated  with  the  greater  tendency  for  striations  to  develop  at  the  higher 
gas  densities. 

A  key  part  of  this  experiment  is  the  measurement  of  the  absolute 
fractional  absorption  from  which  one  can  determine  the  absolute  metastable 
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Figure  9.  Electric  field  to  gas  density  ratio  E/n  for  pulsed  discharges  in 


H2«  The  solid  points  are  our  experimental  data,  while  the  +  and  x 


are  those  of  Breusova  and  of  Guntherschulze,  respectively.  The 


line  is  calculated  from  ambipolar  diffusion  theory  and  theoretical 


ionization  coefficients. 


density.  These  measurements  are  made  using  the  single  mode  dye  laser 
apparatus  and  technique  which  we  developed  previously  for  the  measurement  of 
rare  gas  metastable  densities. Figure  10  shows  a  measured  absorption  line 
profile  for  the  588.8  nm  line  caused  by  transitions  from  the  v  -  1,  N  =  1 
level  of  the  c^nu  state  to  the  v  *  1,  N  =  1  level  of  the  i^n  state. 

Immediately  apparent  is  the  fine  structure  splitting  of  about  6  GHz  for  the 
3  10 

c  iig  state.  Somewhat  bothersome  is  the  apparent  excess  absorption  between 
the  fine  structure  components.  Further  tests  of  such  line  profiles  are 
necessary.  In  order  to  obtain  absolute  c^nu  level  densities  it  is  necessary 
to  know  the  radiative  transition  probability  for  this  line.  At  present  the 
only  data  available  is  the  radiative  lifetime  for  the  upper  state  of  the 

transition. The  radiative  branching  ratio  for  transitions  from  the  i% 

3  3  + 

state  to  the  levels  of  the  c  nu  state  relative  to  the  repulsive  b  Eu  state 

are  unknown.  In  the  absence  of  such  data,  one  use  of  our  results  will  be  the 

determination  of  this  branching  ratio.  Fortunately,  efforts  are  underway 

elsewhere^  to  calculate  radiative  transition  probabilities  for  the  triplet 


system  of  H2 


SECTION  IV 


MODELS  OF  02(a1Ag)  IN  DISCHARGES 

In  this  section  we  summarize  the  current  state  of  our  attempts  to  model 
the  observed  behavior  of  02(a^Ag)  metastables  in  pulsed  02  discharges.  The 
experimental  data  which  we  have  used  for  comparison  is  that  of  Vasileva  £t^ 

I  O 

al.  J  for  a  short  pulsed  discharge  in  Ar-02  mixtures.  These  authors  report 

the  relative  intensities  of  both  the  1.27  pm  emission  from  the  metastables  and 

the  634  nm  emission  resulting  from  collision  of  pairs  of  metastables  as  a 

function  of  02  concentration  in  the  mixture.  More  useful  to  us  is  the 

observation  of  the  relative  time  dependence  of  the  634  nm  dimole  emission. 

The  first  step  in  this  project  was  to  update  our  cross  section  sets  for 

electrons  in  02  and  in  Ar.  This  update  resulted  in  a  set  of  cross  sections 

for  C>2  which  are  much  more  consistent  with  the  electron  beam  data  at  electron 

energies  above  about  20  eV  than  are  the  cross  sections  assembled  in  1976  and 

circulated  in  JILA  Information  Center  Report  No.  28.  The  changes  resulted  in 

negligible  change  in  the  low  energy  discrepancies  pointed  out  in  our  1976 
14 

work.  The  update  of  the  Ar  cross  sections  at  energies  near  excitation 
threshold  was  based  primarily  on  the  recently  published  swarm  experiment  and 
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analysis  of  Tachibana.  A  higher  energies  we  made  use  of  theory  and  electron 
beam  experiments.  As  measured  by  Tachibana,  the  cross  sections  result  in  much 
lower  excitation  rates  for  Ar  raetastables  than  do  previous  cross  section  sets. 

The  second  step  was  the  acquisition  and  adaptation  to  our  problem  of  a 
chemical  kinetics  code  from  the  Chemical  Kinetics  Division  of  NBS.^  Next  a 
set  of  rate  coefficients  for  the  electron  excitation  and  excited  state 
reactions  of  interest  was  assembled.  The  electron  excitation  and  ionization 


rate  coefficients  for  20%  O2  in  Ar  are  shown  in  Fig.  11,  while  the  excited 
state  reaction  rate  coefficients  are  listed  in  Table  III. 


Figure  12  shows  comparisons  of  the  relative  experimental  634  nm  emission 

data  of  Vasileva  et  al.  (solid  curve)  with  our  calculations  of  the  square  of 

the  02(a*Ag)  metastable  density  for  several  different  assumptions  as  to  the 

pulsed  discharge  conditions.  The  dashed  curves  from  our  model  are  labeled 

with  atomic  oxygen  density  and  the  (^(a^Ag)  density  at  the  beginning  of  the 

22  21 

afterglow.  We  note  that  dashed  curve  labeled  1.2  x  10  ,  2.6  x  10  ,  which 

are  the  0  and  (^(a^Ag)  densities  calculated  with  our  electron  excitation  rate 
coefficients,  shows  no  maximum  such  as  found  experimentally.  The  other  dashed 
curves  show  the  afterglow  calculations  when  it  is  assumed  that  no  C^Ca^Ag) 
metastables  are  present  at  the  end  of  the  discharge  and  the  atomic  0  density 


is  varied.  None  of  the  initial  conditions  shown  reproduce  the  time  behavior 
of  the  observed  634  nm  emission.  It  would  appear  that  much  larger  quenching 
rate  coefficients  for  the  O^Ca^A,,)  metastables  or  much  larger  net 
recombination  rates  for  0  atoms  are  required.  Perhaps  the  temperature  rise  of 
about  a  factor  of  two  during  the  afterglow  is  important. 

An  additional  problem  with  our  model  is  that  the  efficiency  of  conversion 
of  0  atoms  into  a*Ag  molecules  is  low,  e.g.,  about  2%,  for  conditions 
corresponding  to  the  discharge  energy  input  in  the  experiments  of  Vasileva  et 
al .  On  the  other  hand,  if  one  assumes  that  all  of  the  input  discharge  energy 


is  used  to  dissociate  O2,  the  calculated  0  atom  density  is  equal  to  the  a^A 


density  which  these  authors  claim  to  be  present  in  the  afterglow  of  their 


discharge.  Some  critical  reactions  in  determining  the  efficiency  of 
metastable  formation  are  the  three  body  formation  of  O3  involving  two  O2 
molecules  and  the  destruction^  of  O3  by  0.  If  both  of  these  reactions  lead 
to  the  formation  of  a*Ag  metastables,  one  begins  to  approach  the  efficiency 
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Figure  11.  Electron  dissociation,  ionization,  and  metastable  excitation 
coefficients  for  33%  O2  in  Ar  vs  E/n  as  calculated  from  our 
revised  electron  cross  section  sets. 
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Reactions  in  afterglow  of  02-Ar  discharge  at  300  K 


TIME  (ms) 

Comparison  of  measured  relative  634  nm  emission  data  with  the 
calculated  square  of  the  C^a^A^)  density  for  a  33!?  0 ^  «  6 7Z  Ar 
mixture.  The  solid  curve  is  taken  from  the  experimental  data  of 
Vasileva  et  al.,  while  the  dashed  curves  are  from  the  model 
discussed  in  the  text. 


required.  Both  of  the  reactions  are  energetically  al Lowed,  but  the  former 
would  have  very  little  excess  energy. 

We  conclude  that  our  model  Is  in  serious  disagreement  with  the  results  of 
Vasileva  et  al.  and  that  the  required  modifications  are  not  obvious.  In  view 
of  the  importance  of  competition  between  0-j  formation  and  a*A g  formation,  it 
would  be  very  desirable  to  have  simultaneous  measurements  of  these  densities. 
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SECTION  V 

CONCLUSIONS 

The  results  described  in  this  report  show  that  all  of  the  metastable 
levels  of  the  c%u  state  of  are  rapidly  destroyed  in  collisions  with 
under  electrical  discharge  conditions.  As  a  result  the  populations  of  the 
metastable  levels  are  expected  to  be  only  marginally  larger  than  those  of 
nearby  radiating  levels,  rather  than  much  larger  as  suggested  in  some  previous 
models  of  these  discharges.  Furthermore,  our  preliminary  measurements  of 
electric  field  strengths  in  the  positive  column  of  pulsed  H2  discharges  show 
that  when  the  energy  input  per  molecule  is  small  the  E/n  values  required  to 
maintain  the  steady  state  discharge  are  very  close  to  values  predicted  using 
our  recently  recommended  electron  collision  cross  section  set.  The  apparent 
success  of  our  techniques  for  electric  field  and  metastable  absorption 
measurements  at  low  energy  input  densities  suggests  future  measurements  at 
higher  energy  densities  in  an  effort  to  define  the  role  of  vibrational ly 
excited  Wj  in  these  discharges. 

The  analysis  of  the  behavior  of  the  (^(a^Ag)  metastables  in  C^-Ar  has 
shown  that  there  is  a  serious  deficiency  in  our  model.  We  are  unable  to 
simulate  the  maxima  in  the  634  nm  emission  produced  by  (^(a^g)  metastables  as 
observed  by  Vasileva  et  al .  Until  this  problem  is  solved  there  does  not  seem 
much  point  in  pursuing  this  modeling. 

In  view  of  the  success  of  the  laser  diagnostic  techniques  in  the 

measurement  of  the  collisional  destruction  rate  coefficients  for  the  ^(c  nu) 

metastable,  we  recommend  that  these  techniques  be  applied  to  measurements  of 

the  collisional  quenching  and  electron  excitation  of  the  higher  metastable 

states  of  N9,  i.e.,  the  a"  ^E  +  and  a'  ^E  “  states.  These  states  have  been 
1  g  u 
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the  subject  of  much  discussion  In  recent  years  because  of  their  possible 
Importance  in  multi-step  ionization  processes  in  N2  discharges  and  shock 
waves. 
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